Generation of Fourier transform limited heralded single photons 
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In this paper we study the spectral (temporal) properties of heralded single photon wavepackets, 
triggered by the detection of an idler photon in the process of parametric downconversion. The 
generated single photons are studied within the framework of the chronocyclic Wigner function, 
from which the single photon spectral width and temporal duration can be computed. We derive 
specific conditions on the two-photon joint spectral amplitude which result in both pure and Fourier- 
transform limited heralded single photons. Likewise, we present specific source geometries which 
lead to the fulfilment of these conditions and show that one of these geometries leads, for a given 
pump bandwidth, to the temporally shortest possible heralded single photon wavepackets. 
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I. INTRODUCTION 

Single photon wavepackets are a crucial ingredient 
for a number of quantum-enhanced technologies such 
as quantum computation with linear optics [1] and se- 
cure quantum key distribution [2 . Recently it has been 
shown that single photons can be efficiently "heralded" , 
whereby photon pairs are generated by the process of 
spontaneous parametric downconversion (PDC) and sub- 
sequently one of the photons is detected, thus condition- 
ally preparing a single photon wavepacket in the conju- 
gate mode [21 m [5l |6j [7]. While great progress has been 
made in the development of so-called on demand single 
photon sources (e.g. based on quantum dots in micro- 
cavities [S], single atoms coupled to optical cavities [S] and 
nitrogen vacancies in diamondflO;), PDC can lead to 
practical, room-temperature sources of heralded single 
photons. A crucial property of PDC-based heralded sin- 
gle photon sources is that the spatial and spectral emis- 
sion properties can be controlled to a large degree by 
appropriate choices of the source geometry, including 
dispersion characteristics of the nonlinear crystal, and 
the pump characteristics. It is indeed remarkable that 
PDC can be configured to produce photon pairs rang- 
ing from those which are factorable to those which are 
highly entangled in continuous degrees of freedom such 
as frequency and transverse wavevector. In this paper we 
exploit this flexibility to derive conditions for the genera- 
tion of both pure and Fourier transform limited heralded 
single photons. Likewise, we present specific source ge- 
ometries which lead to the fulfilment of such conditions. 
Thus, the techniques presented here can be used to de- 
sign PDC photon pair sources leading to heralded single 
photons described by FT-limited, ultrashort wavepack- 
ets. 

Heralded single photon sources rely on the quality of 
the photon-number correlations between the signal and 
idler (used as trigger) PDC modes. It is experimentally 
straightforward to determine such quality by measuring 
the ratio of the coincidence detection rate across the sig- 



nal and trigger PDC modes to the signal singles count 
rate. There are a number of experimental factors which 
can reduce this ratio from its ideal unity value. Higher- 
order PDC contributions whereby more than one photon 
pair can be emitted at a given time implies that a trigger 
detection event (from a non-photon number resolving de- 
tector) can erroneously indicate the presence of a single 
signal photon, while in reality two, or more, photons are 
present. Likewise, the presence of optical losses in the 
signal mode implies that a trigger detection event can 
erroneously indicate the presence of a single signal pho- 
ton, while in reality a vacuum state exists. The presence 
of background photons, for which a trigger is detected 
without the existence of a signal photon evidently also 
reduces the heralding efficiency. In Ref. @] we have pro- 
posed and implemented an experimental criterion which 
assesses source performance, specifically suited to her- 
alded single photon sources, including the effects of op- 
tical losses, higher photon number contributions, back- 
ground photons and the binary behavior of realistic sin- 
gle photon detectors. 

Apart from the quality of photon number correlations 
between the two PDC modes, it is important to consider 
how entanglement in continuous degrees of freedom (e.g. 
spectral and transverse wave-vector) translates into the 
properties of the heralded single photon. Specifically, 
as was shown in Ref. [13], correlations in any degree of 
freedom in the photon pairs results in a departure from 
purity in the heralded single photons. Such impurity is 
undesirable as it implies that multiple single photons pro- 
duced by more than one source will not interfere. In this 
paper we further analyze the spectral (temporal) proper- 
ties of heralded single photon wavepackets, showing for 
a restricted class of states that factorable photon pairs 
lead to both pure and Fourier transform limited heralded 
single photons. We show in addition that in a specific 
source geometry the heralded single photons can attain 
their shortest possible temporal duration, which corre- 
sponds to the pump pulse duration. 
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II. DESCRIPTION OF HERALDED SINGLE 
PHOTONS 

The heralded single photon in the signal mode can be 
represented by its density operator p^: 



Tr,(/5nO, 



(1) 



where Tr^ represents a partial trace over the trigger mode, 
in terms of the density operator of the two photon state 
p and of the measurement operator IIj for the trigger 
(idler) mode: 



Hi = / dujg{uj)\uj)i{uj\i, 



(2) 



which has been expressed in terms of the frequency- 
dependent trigger detection efficiency g{oj). 




FIG. 1: (color online) Heralded single photon source, based 
on pulsed parametric downconversion in a x'^' crystal. 

We will consider a heralded single photon source based 
on the process of parametric downconversion (PDC) 
pumped by an ultrashort pulse, shown schematically in 
Fig. [l] The two-photon state for the PDC process as- 
suming specific directions of propagation (i.e. disregard- 
ing the transverse wavevector degree of freedom) can be 
expressed as: 



I*) 



duJsdljJif{uJs,UJi)\LOs)s\uJi)i 



(3) 



where |w^) = {oj ^)\va,c) ^ with p — s,i in terms of the 
joint spectral amplitude (JSA): 

/(ws,t^i) ^ ^{uJs,UJi)a{uJs + i^i)F{uJs,uJi), (4) 

which is in turn given by the product of the phasematch- 
ing function (PMF) ^(cj^, w^), the pump envelope func- 
tion (PEF) a{uJs + ^j) and a filter function F{ujs,uji). 
While the PMF contains information about the optical 
properties of the nonlinear crystal, the PEF character- 
izes the pump field utilized. The function F(u!s,Wi) de- 
scribes the action of identical spectral filters (labelled ap 



in the figure) acting on both generated photons . The 
JSA is assumed to be normalized so that the integral over 
both frequency arguments of its absolute value squared 
yields unity. Substituting the two-photon density oper- 
ator p — |^')(^'| into Eq. [l] we arrive at the following 
expression for the density operator which results for the 
heralded single photon in the signal mode: 



Ps^ du!ag{uJa)\(l){uja)){(p{uJo)\. 



(5) 



where 



\H^o)) = / duJif{uji,uJo)\uJi) 



(6) 



Let us note than, in contrast to the two-photon density 
operator p = |^')(^'|, in general the single photon density 
operator ps represents a mixed state. In order to gain 
physical insight, let us consider the case where spectral 
filtering is apphed to the trigger mode [such a filter is 
labelled ag in Fig. [l]. In the limit of strong filtering 
we may write the trigger detection efficiency as a delta 
function g{Lo) — > 6{ll! — LUp)- Under these circumstances, 
the resulting single photon density operator is given by: 



Ps = |(/)(wf))(<?!'(wf)|, 



(7) 



which evidently represents a pure state pj]. Thus, an av- 
enue towards a pure heralded single photon source is to 
implement strong spectral filtering on the trigger mode. 
However, this is achieved at the cost of a reduction in 
source brightness as well as a reduction in the bandwidth 
of the prepared photon through non-local state projec- 
tion. 

Is it possible to obtain a pure state without resorting to 
strong spectral filtering? Let us consider the case where 
the joint spectral amplitude f{ujs,u}i) is factorable i.e. 
where functions p{uj) and q{ijj) exist such that f{uJs, uji) — 
p{ujs)q{ii'i)- It is straightforward to show from Eqns. 5] 
and |6] that under these circumstances the heralded single 
photon is indeed pure. In order to gain understanding 
of the conditions for single photon purity, let us perform 
a Schmidt decomposition of the two photon state. The 
JSA can now be expressed as: 



(8) 



/(uJsi^^i) = ^ VX^"m(^s)Wm('^i), 



where Um(w) and Vm{i^) represent the Schmidt functions 
and Am represents the Schmidt eigenvalues |12j. The re- 
sulting two photon density operator ps can be expressed 
as: 



= V KiK J duJog{uJo)Vm{uJo)Vn{^o)\<Pm){4'7^ 



(9) 
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in terms of the Schmidt single photon wavepackets: 



\(f>n 



dLUUm{uj)\Lu) 



(10) 



It is evident from Eq. [9] that in general the heralded 
single photon is described by a mixed state, given as an 
incoherent sum of the Schmidt wavepackets. This expres- 
sion simplifies considerably in the case where the trigger 
detection efficiency is not frequency dependent; using or- 
thogonality of the Schmidt functions, Eq. [9]then reduces 
to: 



5: A. 



J0m)((/' 



(11) 



From an analysis of Eqns. [9] and [TT] it becomes clear 
that if the two photon state contains a single pair of 
Schmidt modes, i.e. if there is a single term in the sum in 
Eq. [8] the resulting heralded single photon state is pure. 
Thus, by engineering the state at the source so as to 
guarantee factorability (for which the state is expressed 
in terms of a single Schmidt mode pair) it becomes possi- 
ble to achieve pure single photon conditional generation 
without resorting to spectral filtering. Conversely, for a 
progressively larger number of active Schmidt mode pairs 
in the photon pair, the attainable single photon purity 
decreases. Indeed, as was shown in Ref. [T5], the single 
photon purity p — Tr(/5g) is given by p = 1/K where 
the Schmidt number K is a, measure of the number of 
active mode pairs. In this case, a pure heralded single 
photon may still be obtained by spectral filtering at the 
cost of a reduction in both single photon bandwidth and 
signal level. Thus, factorable photon pair generation rep- 
resents a crucial enabling technology for future progress 
in experimental quantum information processing. 



III. CHRONOCYCLIC WIGNER FUNCTION 

It is convenient to study the temporal and spectral 
properties of heralded single photons within the frame- 
work of the chronocyclic Wigner function (CWF) formal- 
ism. The latter represents a useful tool for the description 
of classical ultrashort pulses [HJ [TS] and it has also been 
used to represent the PDC two- photon state [TB^. In this 
paper we consider the CWF for the heralded single pho- 
ton, which can be expressed in terms of the single photon 
density operator ps as: 



Ws{co,t) 



1 

2^ 



dLj'{uj + uj'/2\ps\LJ-Lj'/2)e''^'*. (12) 



where subscripts s in the frequency bra and ket have 
been omitted for clarity. Such a function fully charac- 
terizes the spectral and temporal properties of the single 
photons; indeed, the marginal distribution resulting from 
integrating the CWF over frequency or time results in the 
temporal or spectral, respectively, single photon intensity 



profile. In terms of the joint spectral amplitude f{ujs,uii), 
the chronocyclic Wigner function for the heralded single 
photon can be expressed as: 

Ws{uj,t) = ^ f duJogiuJo) I dJ f (uj +^,ujq \ x 



(13) 



Let us note that the chronocyclic Wigner function, 
as expressed in Eq. |13| is a generalization of the 
Wigner- Ville function [T3] for a spectral amplitude func- 
tion /(tJi, 0)2) which depends on two rather than one vari- 
able; the additional integration over the second frequency 
argument results from tracing over one of the two pho- 
tons in the heralding process. 

The spectral intensity profile of the heralded single 
photons is given by the integral over the time variable 
of the chronocyclic Wigner function: 

IU^)= [ dtWs{uj,t)^ /da;og(t^o)|/(^,^o)P, (14) 



while the temporal intensity profile is given by the in- 
tegral over the frequency variable of the chronocyclic 
Wigner function: 



It{t) 



dujWs{uj,t) 



27r 



duJog{uQ) 



dujf{LLi,ujo)e^' 



(15) 



In terms of the Schimdt functions, we may express the 
CWF as: 

Ws{uJ,t) = ^^^V^mK f duJog{^o)Vm{i^o)K{^o) 

m rn 

X J dio'u^niu + ^)<(^ - ^)e-'*. (16) 

If the trigger detection efficiency is not frequency de- 
pendent, using orthogonality of the Schmidt functions, 
Eq. [T6| reduces to: 

Ws{Lo,t) = ^E^™ / y)e^"'tl7) 

In this case, the spectral and temporal intensity pro- 
files are given by: 



m 

It{t) = 2^^\ra\Um{t)[' 



(18) 
(19) 



where Un{t) represents the Fourier transform of u„(a;). In 
order to carry out more explicit calculations, we will as- 
sume PDC photon pairs within a restricted, though quite 



4 



general, class of states. We will consider a frequency- 
degenerate coUinear PDC source (centered at frequency 
oJc) based on type-II phasematching, pumped by a short 
pulse (centered at frequency 2wc) and will disregard 
group velocity dispersion (as well as higher-order disper- 
sion) terms. The phasematching function for such an 
interaction can be expressed as follows: 

(/){ujs,uji) — smc[LAk{ujs,LOi)/2]exp[iLAk{uJs,uJi)/2], 

with the phasemismatch expressed as a Taylor series up 
to first order: 

LAk{Us,UJi) K LAk^-^^ +TsiuJs-UJc) + Ti{uj,-UJc), (21) 

in terms of the PDC central frequency lUc, the constant 
term of the Taylor series AA:'-"-' (which we assume to van- 
ish) and the longitudinal temporal walkofF terms between 
the pump pulse and each of the signal and trigger pho- 
tons: 



Ts = L{kp{2ujc) 



(22) 
(23) 



Let us note that by retaining only constant and group 
velocity terms, we ignore temporal broadening effects re- 
sulting from quadratic and higher order dispersion terms. 
While evidently a complete analysis should include these 
terms, for the case analyzed in detail here which corre- 
sponds to coUinear, degenerate type-II PDC where the 
bandwidth of the generated light is somewhat restricted 
(as compared, for example, to coUinear, degenerate type-I 
PDC), the group velocity approximation is justified. We 
model the pump envelope function as a Gaussian func- 
tion with width <t: 



a{u}s + oji) = exp 



(24) 



In order to obtain an analytic expression for the CWF 
in closed form, we approximate the sine function appear- 
ing in the phasematching function (see Eq. 20 1 as a Gaus- 
sian function [sinc(a;) ~ exp(— 72;^) with 7 w 0.193]. As 
we will discuss below, comparisons of the CWF calcu- 
lated numerically without resorting to approximations 
and the analytic CWF suggest that the Gaussian ap- 
proximation represents a valid approximation in that it 
yields a useful indication of the single photon bandwidth 
and temporal duration. Likewise, we assume a Gaussian 
shape (with central frequency ujgQ and bandwidth Ug) for 
the trigger detection efficiency g{uj) [see Eq. |5]: 



g{uj) = exp 



(25) 



and a Gaussian shape for the two-photon filter function 
FiuJs-i'^i) with width up'- 



F{ujs,uji) = exp 



(26) 



Carrying out the integrals we obtain a CWF of the 
following form: 



1 



nAujAt 



exp 



AJ^ 



X exp 



Ai2 



(27) 



in terms of the temporal and spectral widths At and Alo, 



At = T,. 



Auj 



T ■ 



y ss ii si 

temporal and spectral shift terms T and il: 
T = Ts/2, 

1 T2 

g ss ii si 



(28) 
(29) 

(30) 
(31) 



where we have defined the quantities T^^, and Tgi as 
follows: 



rpZ 
SS 
rji2, 

^ ii 

rp2 



= 2/4 + 2/C72 + ^r2/2, (32) 
- l/r72 + 2/4 -h 2/^2 +7rf/2, (33) 
= 2/^2 + ^r,r,/2. (34) 



It is straightforward to prove that for a two photon 
state for which its corresponding joint spectral ampli- 
tude can be expressed fully in terms of Gaussian func- 
tions, the spectral width Aw and temporal duration At 
must satisfy the inequality At Auj > 1; here Aw and At 
represent half widths at e~^. Note that this is a par- 
ticular instance of the uncertainty relation StScu > 1/2 
which is fulfilled for a broad class of functions (where 
6 refers to the standard deviation). However, the two 
photon state is given in terms of a sine function, which 
does not have a well defined variance and therefore we 
cannot apply this result. Our approach is rather to use 
the Gaussian approximation in cases where it yields es- 
sentially the same chronocyclic Wigner function as with- 
out recourse to such an approximation, and employ the 
uncertainty relation associated with Gaussian functions. 
In order to study possible regimes resulting in Fourier 
transform-limited conditional single photon generation, 
we compute the time-bandwidth product TB = AcuAt: 



TB = A<Aw = (1 - T^JT^^T^ 



(35) 



The quantity TB assumes its minimum value (TB = 1) 
for Fourier transform limited single photons. Such pho- 
tons are characterized by the important property that 
their temporal duration At is the shortest possible com- 
patible with a given spectral width Aw. For single pho- 
tons which are not Fourier transform limited, the time- 
bandwidth product assumes a value greater than unity; 
indeed, the specific value of TB provides a convenient 
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measure of the departure from the Fourier transform 
limit. Clearly, the greater the numerical value of TB, 
the longer the time duration At = TB/Ao; for a given 
fixed spectral width. 

Besides the temporal duration of the heralded single 
photon, another important characteristic time to con- 
sider is the correlation time: the time of emission differ- 
ence between the signal and idler photons. Such a time 
is important because it defines the uncertainty in the ex- 
pected detection time of the heralded photon relative to 
the detection of the trigger photon. In order to calcu- 
late the correlation time Tc, we first compute the joint 
temporal intensity where f{ts,ti) denotes the 

two-dimensional Fourier transform of the joint spectral 
amplitude f{ujg,uji). The joint temporal intensity rep- 
resents the probability density function for generating a 
photon pair at specific times tg and ti. As a second step, 
we may express the joint temporal intensity in terms of 
the variables t+ = tg + U and i_ = ts — ti. From this, 
we can calculate the marginal probability distribution 
S-(t-) which yields the probability density function for 
generating a photon pair with a certain time of emission 
difference i_ . For the approximations already introduced 
in this paper, we may express this t_ marginal distribu- 
tion as S-{t^) — Tr^^/^T^^exp{—t^/T^) in terms of the 
correlation time r„: 



+ -faliTj - T,y 
V2aF 



(36) 



which in the absence of spectral filtering (ap — > oo) re- 
duces to: 



where: 



(37) 



(38) 



It is instructive to compare Eq. 28 (together with 32 ) 



for the heralded single photon time duration with Eq. |37| 
(together with Eq. 38 1 for the correlation time. It be- 



comes apparent that while the former is defined by a 
group velocity mismatch term between the pump and 
the signal photon, the latter is defined by a group veloc- 
ity mismatch term t_ between the trigger and the signal 
photons. It should also be stressed that the correlation 
time does not depend on the pump bandwidth [T7]. 



IV. FT-LIMITED SINGLE PHOTONS VIA 
GROUP VELOCITY MATCHING 

In order to make the discussion more specific we will 
now consider the concepts developed so far for particular 
experimental situations. Let us first consider the case 
of PDC generation in the limit of vanishing pump band- 
width (i.e. a 0). It is clear from Eqns. [28]and 32 that. 



the single photon time duration becomes dependent only 
on the pump bandwidth: At = V2/a. As might be 
expected this time duration tends to infinity in the ide- 
alized CW-pump case. Likewise, the spectral width in 
the cr — > limit can be computed from Eq. 29 and yields 
Auj = ■\/2/7/|r_ I [which in fact corresponds to the un- 
filtered reciprocal correlation time (see Eq. 37 1]. Thus, 



the time bandwidth product in the limit of a very short 
pump bandwidth TB = 2/{y/ja\T-\) diverges. We con- 
clude therefore that CW-pumped PDC light is ill-suited 
for generation of Fourier transform limited heralded sin- 
gle photons. 

Let us briefly consider degenerate, coUinear type-I 
PDC as a candidate for ultrashort heralded single pho- 
ton generation. Within the first order treatment pre- 
sented in this paper (ignoring group velocity dispersion 
and higher order dispersive terms in the phasemismatch 
Ak{LL!s,LL!i)), degenerate coUinear type-I photon pairs are 
characterized by = or equivalently by t_ =0. It has 
been shown [TB] that the orientation of the phasematch- 
ing function in {u!s,Wi} space is given by — arctan(Ts/ri), 
so that in the case of degenerate coUinear type-I PDC, 
this orientation is fixed at —45°. Furthermore, within 
the linear approximation used, this function has an infi- 
nite extent along the direction ujs — uji and shows perfect 
overlap with the pump envelope function which is char- 
acterized by the same orientation: i.e. all frequency pairs 
within the pump envelope function are correctly phase- 
matched. Thus, within this approximation the spectral 
width of each of the two single photon wavepackets is 
infinite, while the temporal duration has a finite value, 
determined by the pump bandwidth and crystal length. 
Let us stress that the inclusion of GVD terms is nec- 
essary in the degenerate, coUinear type-I case in order 
to obtain a physically meaningful result (while generally 
speaking the linear approximation is justified for degen- 
erate, coUinear type-II PDC), leading to a very large 
(but finite) spectral bandwidth for each photon. The 
very large spectral bandwidth (which becomes extraor- 
dinarily large when the PDC photons are near the zero 
group velocity dispersion wavelength[19 ), however, un- 
fortunantely does not result in ultrashort single photon 
wavepackets (though it does imply a very short correla- 
tion time, see discussion above and Eqns. 36 and|37|. In 
Fig. [2] we present the chronocyclic Wigncr function (cal- 
culated numerically from Eq. 13 1 for a specific coUinear, 



when the pump bandwidth is much smaller than the re- 
ciprocal group velocity mismatch term (i.e. a <ti 1/ts), 



degenerate type-I example. For a 5mm long BBO crys- 
tal (with a cut angle of 29.2°), pumped by an ultrashort 
pulse train with 5nm pump bandwidth producing pho- 
ton pairs centered at 800nm we obtain At = 288.01fs 
and Aui = 61.4THz, corresponding to AtAu « 17.7[T5] 
which suggests a large departure from the transform 
limit. Apart from the above considerations, let us note 
that degenerate, coUinear type-I signal and idler photon 
pairs can at best be split probabilistically leading to a 
maximum heralding efficiency of 50%. Thus, we like- 
wise conclude that degenerate, coUinear type-I PDC is 
ill-suited for generation of Fourier transform limited her- 
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aided single photons. 
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FIG. 2: Chronocyclic Wigner function (CWF) obtained nu- 
merically for a BBO crystal of 5mm length, phasematched for 
coUinear, degenerate type-I PDC centered at SOOnm, pumped 
by an ultrashort pulse train with a 5nm bandwidth. The re- 
sulting single photon time duration At — 288.01fs and single 
photon bandwidth Au; = 61.4THz, with AtAuj ^ 17.7, sug- 
gest a large departure from the transform limit. 

Let us turn our attention to the process of coUinear, 
degenerate type-II PDC pumped by ultrashort pulses. In 
the propagation of short optical pulses through dispersive 
media, group velocity dispersion and higher dispersive 
terms result in temporal broadening, while group veloc- 
ity terms result in propagation without distortion. The 
situation is somewhat different for the generation of PDC 
light by short pump pulses, where group velocity terms 
in the nonlinear crystal can result in photon pairs where 
the constituent single photon wavepackets are temporally 
much longer than the pump pulse. We may write down 
the heralded single photon duration from Eq.|28]in terms 
of the pump pulse duration Tp — V2/a as: 



At = 




(39) 



It is evident from Eq. [39] that the shortest possible 
heralded single photon temporal duration is the pump 
pulse temporal duration itself. Likewise, we conclude 



from Eq. 39 that this shortest temporal duration is ob- 
tained in the limit of no spectral filtering ((jp ^ oo) and 
of negligible group delay as compared to the pump pulse 
duration i.e. Ts <^ Tp. It becomes apparent that for 
an ultrashort pump pulse, the attainable temporal dura- 
tion is limited by the group velocity mismatch term r^. 
Two possibilities exist for making r^/r^ ^ 1. The first 
involves using a very short crystal, concretely one obey- 
ing L <C Tp/{k'p — kg); indeed for a short enough crys- 
tal the group delay due to group velocity mismatch 
can become insignificant. This involves the need for pro- 
gressively shorter crystals as the pump pulse duration is 



decreased. It clearly does not represent a viable solu- 
tion as it implies a prohibitive reduction in PDC flux. 
The second possibility is to employ a PDC interaction 
obeying group velocity matching between the pump and 
the PDC light. Indeed, it is clear from Eq. 39 that for 



= (which implies k'p — k'^) we can obtain, with no 
spectral filtering, the shortest possible time duration Tp 
irrespective of the crystal length used. 

Let us stress that in our analysis in this paper, we 
have ignored group velocity dispersion terms (as well as 
higher order dispersive terms). For classical pulse prop- 
agation, under these circumstances we would expect no 
temporal broadening. However, Eq. [35] for the heralded 
single photon time-bandwidth product does indeed per- 
mit values greater than unity, even for an unchirped 
pump. The physical mechanism by which the single pho- 
ton wavepackets which constitute a PDC photon pair suf- 
fer temporal elongation is related to group velocity mis- 
match between the pump pulse and the generated PDC 
light. Indeed, in the long crystal limit the heralded single 
photon duration is determined purely by the group veloc- 
ity mismatch term: Atg — > y/jTs/V2. In order to under- 
stand this physically, let us consider a photon pair cre- 
ated on the second face of the crystal and secondly a pho- 
ton pair created on the first face. While in the latter case 
the amplitudes for photon pair generation emerge from 
the crystal temporally overlapped with the pump pulse, 
in the former due to group velocity mismatch, these am- 
plitudes emerge ahead of the pump pulse by a time Tg. 
Because the wavepackets corresponding to each of the 
two photons are the result of adding the amplitudes from 
all locations within the crystal, the result is a temporally 
stretched wavepacket (as compared to the pump pulse). 
As a numerical example [to be discussed further below; 
see Fig. |3jC)], let us consider a 5mm long BBO crystal 
pumped by a pulse train centered at 400nm with 5nm 
bandwidth. The expected duration of the wavepackets 



computed from Eq. 28 is 127. 2fs (for a pump pulse with 



28.3fs duration), resulting in a time-bandwidth product 
of TB = 6.4. As has already been discussed, for small 
pump bandwidths together with type-1 operation, the 
time bandwidth product can take much larger values; in- 
deed it can be shown that values as high as 10^ for the 
Schmidt number K and TB are possible [IQj. 

How can we achieve Fourier transform limited condi- 
tional single photon generation? An analysis of Eq. [35] 
reveals that TB = 1 if the correlation coefficient 
Tg^/ {TssTii) vanishes. One way in which this can occur 
is if strong filtering is applied to both photons {ap ^ Q) 
or if strong filtering is applied to the trigger photon only 
{(jg — > 0). In both of these cases, it can be readily shown 
that in the limit of strong filtering, the Fourier transform 
limit is attained (i.e. TB — > 1), evidently at the cost of 
a prohibitive reduction in flux. In addition, it should 
be stressed that while the FT limit may be attained via 
filtering, the bandwidth is thereby reduced thus preclud- 
ing ultrashort single-photon wavepackets. Alternatively, 
the correlation coefficient T'^^/ {TssTu) vanishes if Tsi = 
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which in fact represents the condition derived in Ref. (2D] 
for spectral factorizability in photon pair states. Thus, 
we reach the important conclusion that factorable two 
photon states [131 [201 HIl US] lead to Fourier transform 
limited heralded single photons. Let us note that the 
latter route involves no filtering, and therefore can lead 
to bright, broadband sources of FT- limited single pho- 
tons. 



An analysis of the CWF [see Eq. 12 reveals that the 
single photon spectrum is in general shifted from the 
central PDC frequency uJc by a frequency [given by 
Eq. |31| which depends linearly on the frequency detun- 
ing of the trigger filter with respect to the central PDC 
wavelength. Therefore, as the trigger filter passband is 
shifted, the resulting heralded single photon spectrum 
also shifts. This is a direct consequence of spectral cor- 
relation in the photon pair state; filtering of one of the 
two photons, non-locally projects the conjugate photon 
into a particular spectral band. Indeed, U, is proportional 
to Tgi , which controls the degree of factorizability in the 
photon pair [201. Thus, for a factorable two photon state 
(for which Tsi — 0) the heralded single photon spectrum 
does not shift due to shifts of the trigger filter passband. 
This suggests an experimental test for the degree of fac- 
torizability (and thus for the degree of departure from the 
FT limit): monitoring the heralded single photon central 
frequency as a function of the trigger filter central pass- 
band frequency. 



The factorizability condition Tsi — (see Eq. 34 1 can 
be satisfied, for example, for the symmetric group veloc- 
ity matching (SGVM) condition [20 1 l22 l |23]. expressed as 

-\- Ti = Q, ov alternatively as: 



2k'{2ojc)-K{uJc)-K{uj,)^{). 



(40) 



Physically, this condition implies that the signal and 
idler photon amplitudes walk away longitudinally from 
the pump in a symmetric fashion: i.e. one outpaces the 
pump pulse by the same group delay as its conjugate lags 
the pump pulse. The resulting single photon temporal 
duration can be computed from Eq. [28] together with the 
conditions Ts + ti = and Tgi — (see Eq. 34 1 , obtaining 
the result 



M = \/2r„ 



Thus, in the symmetric group velocity matching sce- 
nario, the single photon time duration assumes, within 
a factor of \/2, its smallest possible value. It may be 
similarly shown that the spectral width is in this case 
simply the reciprocal of the temporal duration so that 
the Fourier transform limit is attained without recourse 
to filtering. 

Let us study an alternative regime for Fourier trans- 
form limited single photon generation, where the pump 
group velocity is matched to the signal photon group ve- 
locity: i.e. fcp(2wc) = fcs(wc) or = 0, a regime which we 



refer to as asymmetric group velocity matching (AGVM) . 
In this case, it may be readily shown that TB — > 1 if 
CTTi ^ 1. The latter condition may be re-expressed as: 



(42) 



Thus, by imposing the asymmetric group velocity 
matching condition together with a large crystal length- 
pump bandwidth product, it becomes possible to gen- 
erate Fourier transform limited heralded single photons. 
AGVM is preferable over SGVM in the sense that for the 
former the heralded single photons attain the shortest 
possible time duration At = Tp (equal to the pump du- 
ration) together with the largest FT-limited bandwidth 
possible Aix) — \/Tp (while for the latter these parameters 
differ from the ideal values by a factor of v^). However, 
it should be stressed that AGVM leads to signal and idler 
photons which are distinguishable from each other; the 
former is not a limitation for heralded single photon gen- 
eration but would be a limitation for experiments where 
the signal and idler photons are made to interfere with 
each other. Indeed, an advantage of the SGVM technique 
is that it combines ultrashort FT-limited heralded single 
photon generation with a symmetric joint spectral ampli- 
tude yielding indistinguishable signal and idler photons 
(except in polarization). 

Another important aspect that should be considered 
when designing a heralded single photon source is the 
correlation time. This time determines the uncertainty in 
the expected time of arrival of the heralded single photon 
with respect to the trigger detection time. As has already 
been discussed, type-I PDC can result in extremely short 
correlation times, due to the larger PDC bandwidths pos- 
sible. However, as has also been discussed, despite the 
short correlation times, type-I leads to states exhibiting a 
strong deviation from the FT limit. Let us compare her- 
alded single photon sources resulting from SGVM and 
AGVM in terms of their correlation times. For SGVM, 
Tg — —Ti and therefore Eq. 



37 



reduces to Tc — y/2j\T, 



while for AGVM Ts = and therefore Eq. [37] reduces 
to Tc — y^j/2\Ti\. When comparing these two results, it 
should be recalled that depends linearly on the crys- 
tal length and that in order to obtain a factorable state 
through AGVM, it is essential to use fairly long crystals 
What this means in practice is that the cor- 



(41) [seeEq.|42 



relation time tends to be much longer for heralded single 
photons produced through AGVM compared to heralded 
single photons produced through SGVM. 

Let us consider specific examples to illustrate the single 
heralded photon duration and correlation time behavior. 
First, let us consider a 2cm KDP crystal cut at 67.7° 
(coUinear, degenerate type-II emission at 830nm) which 
can yield FT-limited single photons through AGVM [13]. 
Fig. [3| A) shows the CWF computed numerically for this 
case, assuming a pump duration Tp = 30.4fs (which cor- 
responds to a 5nm pump bandwidth) . The resulting sin- 
gle photon duration under the Gaussian approximation 
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is At = Tp (where Aui = 1/At), while the correlation 
time is 896. 3fs. Secondly, let us consider a 2.3nim-length 
BBO crystal cut at 28.8° (coUinear, degenerate type-II 
emission at 1514nm) which can yield FT-limited single 
photons through SGVM[20^. Fig. [sJB) shows the CWF 
computed numerically for this case, assuming a pump du- 
ration Tp — 33.7fs (which corresponds to a 15nm pump 
bandwidth). The resulting single photon duration un- 
der the Gaussian approximation is At = \f2Tp (where 
Auj = 1/At), while the correlation time is 67.5fs (note 
this value is independent of the pump bandwidth used) . 
Thus, while AGVM leads to the shortest possible time 
durations, it involves distinguishable signal and idler pho- 
tons and longer correlation times as compared to SGVM. 

In Fig. |3]^C) we show, for comparison, the CWF cor- 
responding to a PDC source which docs not fulfil any 
particular group velocity condition and which therefore 
cannot be used as the basis for an FT-limitcd heralded 
single photon source. Here, we have assumed a 5mm- 
length BBO crystal cut at 42.3° for coUinear, degener- 
ate type-II emission at 800nm. The resulting single pho- 
ton duration obtained numerically is At = 205. 7fs while 
the single photon bandwidth is Aw = 50.1THz, yielding 
AtAio K, 10.3[TS], which suggests a substantial devia- 
tion from the Fourier transform limit. Figs. [3|D)-(F) 
show the contour defined by frequencies and times for 
which the CWF takes a value of times the maxi- 
mum value, calculated numerically (solid line) and an- 
alytically through Eqn. [27] (dashed line) for each of the 
examples in Figs. [3](A)-(C). In general terms, while the 
single photon bandwidth calculated numerically and an- 
alytically agree very well, the analytic result based on 
the Gaussian approximation tends to underestimate the 
temporal width. A computation of the temporal widths 
calculated numerically At„„ and analytically Atan while 
varying source parameters, such as crystal length and 
pump bandwidth (not shown) reveals that At„„/Atan 
ranges in value roughly from 1 to 1.6. This is due to the 
difficulty of approximating top-hat-type functions (which 
result from a sinc-type structure in the spectral domain) 
with Gaussian functions. It should also be pointed out, 
however, that there is a class of PDC two-photon states 
for which the Gaussian approximation leads essentially 
to the same result as predicted without resorting to this 
approximation. This includes those states for which the 
crystal length L or the pump bandwidth a (or both) 
are small. Furthermore, the maximum product Ltr for 
which the Gaussian approximation can be used reliably 
depends on the orientation of the PMF in signal - idler 
frequency space. Thus, for cases where the PMF is ori- 
ented vertically such as the AGVM example above, the 
sine structure is fully contained by the idler single photon 
wavepacket, which is traced over by heralding so that it 
has no effect on the prepared signal photon wavepacket. 
The resulting agreement between the numerical and ana- 
lytic calculations [see Fig.[3]jD)] is essentially perfect; the 
very small disagreement is due to quadratic and higher 
order dispersion terms which are ignored by the analytic 



approach. Likewise, for the SGVM example above, the 
crystal is sufficiently short that the agreement is very 
good though the discrepancy is somewhat more apparent 
than for the AGVM example [see Fig.|3jE)]. Finally, for 
the third example shown, for which no particular GVM 
condition is fulfilled, the disagreement is more apparent. 
In this case while the analytic and numerical calculations 
both predict Auj — 50.1TIIz, the time durations obtained 
differ as follows: Atan = 127. 2fs and At^u — 205. 7fs [see 
Fig-HtF)]. 
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FIG. 3: (color online) Panels (A)-(C) show examples of nu- 
merically calculated chronocyclic Wigner functions (CWF) 
without resorting to the Gaussian approximation and tak- 
ing into account all dispersive orders, centered at a time t=0, 
for a single heralded single photon source based on: (A) a 
KDP crystal, fulfilling asymmetric group velocity matching, 
(B) a BBO crystal, fulfilling symmetric group velocity math- 
cing, and (C) a BBO crystal not fulfilling any special group 
velocity condition. Panels (D)-(F) show, correspondingly to 
(A)-(C), the contour associated to one half of the maximum 
value for the CWF calculated numerically (solid line) and for 
the CWF calculated analytically (dashed line; see Eg. |27[ ). 

We have shown that a factorable two-photon state 
leads to heralded single photons which are both pure and 
FT limited. We now address the question of what the re- 
lationship is, in general, between the Schmidt number 
K (which quantifies the departure from photon pair fac- 
torability) and the time-bandwidth product TB (which 
quantifies the departure for the single photons from the 
FT limit). Under the Gaussian approximation for the 
phascmatching function it is possible to carry out an an- 
alytical Schmidt decomposition of the two photon state 
expressed in Eq. [s] (together with Eq. |4|; such a cal- 
culation was carried out in Ref. [20 . For balanced fil- 
tering for both photons (og co). it can be readily 
shown that the expression for the Schmidt number thus 
obtained is exactly identical to the expression we have 



derived for the time-bandwidth product TB [see Eq. 35 
We therefore complement the results previously derived 
that the single photon purity p is the reciprocal of the 
Schmidt number [13] and in turn equal to the visibility 
V expected in a so called event-ready Hong-Ou-Mandel 
(HOM) interferometer |24j (where two identical heralded 
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single photons, from separate sources, are made to in- 
terfere in a HOM geometry) with the new result that 
the time-bandwidth product is also equal to the Schmidt 
number. In summary, we arrive at the following equation 
linking the degree of spectral entanglement in the pho- 
ton pairs (K) with important quantities relating to the 
resulting heralded single photon: 

K ^TB^l/p^l/V. (43) 

While recent experiments demonstrate nearly ideal sin- 
gle photon heralding efficiencies, the generation of spec- 
trally pure, FT-limited single photon wavepackets repre- 
sents a challenge which must be addressed to make such 
sources useful for quantum information processing appli- 
cations. In this context, the relationships expressed by 
Eq. |43] are important for the design of optimal heralded 
single photon sources. Eq. |43] relates an important prop- 
erty of the photon pairs used (the degree of spectral en- 
tanglement, or Schmidt number K) with crucial proper- 
ties of the resulting heralded single photons. These rela- 
tionships suggest possible experimental procedures to de- 
termine source quality. Specifically, a measurement of the 
visibility in an event ready HOM interferometer would 
yield both the purity and the time-bandwidth product of 
the interfering single photon wavepackets. 

V. CONCLUSIONS 

We have studied within the chronocyclic Wigner func- 
tion (CWF) formalism the temporal (spectral) properties 
of heralded single photon sources based on the process 
of parametric downconversion, specifically in a type-II, 
coUinear, frequency-degenerate regime. We have estab- 
lished, through a Gaussian approximation of the phase- 
matching function, a clear relationship between the prop- 



erties of the photon pairs generated by coUinear, degener- 
ate type-II PDC and those of the resulting heralded single 
photon. Despite the Gaussian approximation used, for an 
important subset of possible states, the CWF predicted 
analytically is essentially identical to that obtained nu- 
merically without recourse to approximations. We have 
shown that spectral entanglement in the photon pairs, as 
quantified by the Schmidt number, translates into spec- 
tral mixedness and a departure from the FT-limit in the 
heralded single photons. Conversely, we have shown that 
factorable two photon states lead to both spectrally pure 
and Fourier transform limited heralded single photons. 
We have furthermore discussed two distinct group veloc- 
ity matching scenarios both of which lead to factorable 
two photon states and therefore to pure and Fourier 
transform limited heralded single photons, without re- 
course to filtering. We have shown that heralded single 
photons obtained via asymmetric group velocity match- 
ing lead to the shortest possible single photon time du- 
ration, which corresponds to the pump pulse duration. 
The former is achieved, however, at the cost of a rela- 
tively large correlation time and of distinguishable signal 
and idler photons. Alternatively, heralded single photons 
based on symmetric group velocity matching lead to a 
single photon duration greater than the ideal value by a 
factor of \/2 while exhibiting signal-idler indistinguisha- 
bility and the lowest possible FT-limited correlation time. 
We believe these results may be important for the fur- 
ther development of practical single photon sources for 
quantum information processing. 
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